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Abstract. In this work we introduce and analyze a linear size-structured pop- 
ulation model with infinite states-at-birth. We model the dynamics of a pop- 
ulation in which individuals have two distinct life-stages: an "active" phase 
when individuals grow, reproduce and die and a second "resting" phase when 
individuals only grow. Transition between these two phases depends on in- 
dividuals' size. First we show that the problem is governed by a positive 
quasicontractive semigroup on the biologically relevant state space. Then we 
investigate, in the framework of the spectral theory of linear operators, the 
asymptotic behavior of solutions of the model. We prove that the associ- 
ated semigroup has, under biologically plausible assumptions, the property of 
asynchronous exponential growth. 
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1. Introduction 



Interest in understanding the dynamics of biological populations is old. Clas- 
sical, ordinary differential equation models assume homogeneity of individuals 
within population classes, and involve equations for total population sizes. How- 
ever, individuals in biological populations differ in their physiological character- 
istics. Vital rates, such as those of birth, death and development, vary amongst 
individuals. Therefore physiologically structured partial differential equation mod- 
els are often more useful to understand the dynamics of biological populations. We 
refer the interested reader to the monographs [T71 [5D1 HT] for basic concepts and 
results in the theory of structured populations. 
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In this paper we study the foUowmg Imear size-structured model 

ui^tis,t) + {'ji{s)ui{s,t))^ -H{s)ui{s,t) + / I3{s,y)ui{y,t)dy 

Jo 

~ Ci{s)ui{s,t) + C2{s)u2{s,t), (1.1) 
U2.tis,t) + (72(5)1X2(5, t))^ = Cl{s)ui{s,t) - C2(s)u2(s,t). (1.2) 

Equations (jl.ip - (|1.2p are equipped with the following boundary and initial condi- 
tions 

7i(0)ui(0,t) =0, ui{s,0)^ul{s), 

72(0)^2(0, t) = 0, u2{s,o) = u;{s). ^ • ' 

Individuals may experience two different stages in their life that we call "active" 
and "resting" . The density of individuals in the active stage of size s G [0, m] 
at time t is denoted by ui(s,t), while the density of individuals in the resting 
stage of size s at time t is denoted by U2{s,t). The maximal size an individual 
may reach is denoted by m. Individuals grow in both classes at the size-dependent 
growth rates 71 and 72, respectively. In the active stage individuals experience 
size-dependent mortality denoted by fx. Further, only individuals in the active 
stage reproduce, this is expressed via the recruitment term on the right hand side 
of equation p.ip . In particular, the function (3{s,y) gives the rate at which an 
individual of size y produces offspring of the size s. The transition between the 
two life-stages is captured by the size-dependent rates ci and C2. We make the 
following assumptions on the model parameters 

Ai, ci, C2 eL^([0,m]), 0</3GC([0,m]2), (3^0, < 71, 72 e ^^([0, m]). 

(1.4) 

Note that our assumptions in (|1.4[) are tailored toward the mathematical analysis 
of model (|l.ip - p.3p . In case of a specific population, one can make additional 
assumptions on the model ingredients, such as P{s,y) ~ whenever s > y, that 
is, individuals can only produce offspring of smaller size. For later use, let 

B=\\j3\\oo and C = max{||ci||oo, ||c2||oo}, 

where ||.||oo stands for the usual norm. In [3j and 0, the authors proposed 
and studied a similar linear age-structured model that describes the dynamics 
of a population which consists of proliferating and quiescent cells. It was shown 
that under some conditions on the support of the transition rates between the 
compartments, the semigroup associated with the linear problem has the property 
of asynchronous exponential growth [31 [7] ■ In contrast, in our model individuals 
are structured by size, individual development depends on size. Another differ- 
ence is that we consider a very general type of recruitment term (see e.g. [3]), 
namely we assume that individuals may have different sizes at birth. We also refer 
the interested reader to and |16] where one-dimensional linear size-structured 
cell population models, with different recruitment terms, were investigated in the 
framework of positive operator theory. 
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2. Existence and positivity of solutions 

Our main objective in this section is to show that our model (|l.ip - p.3p is 
governed by a strongly continuous (Co for short) positive quasicontractivc semi- 
group on the biologically relevant state space X = L^(0, m) x L^(0,m). This will 
imply that problem p.ip - (|1.3p is well-posed and every solution starting with a 
non-negative initial condition remains non-negative. 

We equip the space X with the norm ||x||;\;' — ||a;i||i -I- ||a;2||i, where x = 
{xi,X2) and || • ||i stands for the usual norm. For x ~ {xi, X2), y = (2/1,2/2) S X 
we define the partial ordering x < y if and only if xi{s) < yi{s) and X2{s) < 2/2(5) 
for a.e. s £ [0, m]. Then X is a. Banach lattice. We refer to the book by Arendt 
et al. O part C] for definitions and basic facts about Banach lattices. The dual 
space of X is X* = L°°(0,m) x L°°(0, m) equipped with the norm ||x*||;f. = 
sup{||a;|||oo> ||2;2llc5o}: where x* = {x\,X2) and || ■ ||oo stands for the usual L°° 
norm. Finally, the natural pairing (or semi-inner product) between elements x = 
{xi,X2) G X and x* = {x^^xD G X* is given by 

(x,x*)a'= / {xi{s)xl{s) + X2{s)x*2{s))ds. 
Jo 

Next we cast system (|l.ip - (jl.3p in the form of an abstract Cauchy problem on the 
state space X as follows 



-u={A + B)u, u(0) = uo, 



(2.1) 
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where u = (mi,M2). The unbounded part is given by 

with dense domain 

Dom(^) = {u G W^-^{0,m) x W^-^{0,m) \u{0) = O} . 

Notice that here is the zero vector in M^. The bounded part is given by 

f-ifJ' + ci)ui+ J^" P{-,y)ui{y)dy + C2U2\ ^ 
\ -C2U2 + C1U1 J 

Our aim is to estabhsh that the hnear operator ^ + B is a generator of a positive 
quasicontractive semigroup (see [UITO]). To this end first we recall (see e.g. [HEIIH]) 
some basic concepts and results from the theory of linear operators acting on (non- 
reflexive) Banach spaces. Let £ be a linear operator defined on the real Banach 
lattice y with norm \\.\\y. C is called dissipative if for every A > and y G Dom(£), 

\\{I-XC)y\\y>\\y\\y. 

A function / : 3^ ^ M is called suhlinear if 

f{y + z)<f{y) + f{z). y.z^y 

f{Xy)^\f{y), A>0, yey. 

If in addition f{y) + f{—y) > holds true for y ^ then / is called a half-norm 
on y. The linear operator C is called f- dissipative if 

f{y)<fiy-XCy), A > 0, yeBom{C). 

An operator C is called dispersive, if it is p-dissipative with respect to the canonical 
half-norm 

P{y)^\\y+\\y, (2.2) 

where y^ ~ 2/ VO (and y~ ~ (— y)+). In our setting, the positive part of an element 
X = (xi, X2) of the state space X is defined as follows 

{xi{s),X2{s)) if xi(.s) > 0, X2(s) > 0, 

(2;i(s),0) if .Ti(.s) > 0, a;2(s) < 0, , . 

iO,X2{s)) if xi{s) < 0, X2{s) > 0, • 

(0,0) if .Ti(s) < 0, X2is) < 

Clearly, p defined by (|2.2p is a continuous sublinear functional and its subdiffer- 
ential is given by 

dp(y) = {(/>y G yi such that \\(t)y\\y' < 1, {y,(t>v)x ^\\y'^\\y) , 

where is the positive cone of y* . We also note that it follows from the Hahn- 
Banach Theorem that dp{y) ^ for every y S 3^. In fact, (j)y G dp (y) if and only 
if (l)y{s) = 1 if y{s) > 0, < (l)y{s) < 1 if y{s) = and 0y(s) = if y{s) < 0. 
We recall that a Co semigroup T(t) is called quasicontractive if 

||r(t)|| <e"*, t>0. 
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for some G R, and it is called contractive if a; < 0. Quasicontractive semigroups 
on spaces are of special interest, for example it can be shown (see e.g. [T5] ) 
that every quasicontractive semigroup on an space admits a minimal domi- 
nating positive semigroup (so called modulus semigroup, see e.g. [51 p. 278] for a 
definition), that itself is quasicontractive. Let us first recall the following charac- 
terization theorem from [5l Corollary 7.15] (sec also Theorem 1.2 in Sect. C-II in 
0). 

Theorem 2.1. Let y be a Banach lattice and let C : Dom{C) y be a linear 
operator. The following statements are equivalent. 
(i) C is the generator of a positive contractive semigroup. 

(it) C is densely defined, Rg{XX — C) = y for some A > 0, and C is dispersive. 

We also recall (see e.g. [2]) that £ is dispersive if for every y G Dom(£) there 
exists 4>y G y* with < (j)y, \\4'y\\y* < 1 and {y,(f>y)x — Wu^Wy such that 

{Cy,cby)x<0- 

In fact, we have (jjy G dp (y). Let us finally recall the following generation theorem 
from fH Corollary 5.8]. 

Theorem 2.2. Let T{t) and S{t) be strongly continuous semigroups on X with 
generators A and B, respectively, satisfying 



<Me'^\ for all t > 0, n G N, 

for some constants M > 1 anduo G M. Let D = Dom{A)C]Dom{B) and assume that 
D and (Ao — A — B)D are dense in X for some Aq > OJ. Then C — A + B generates 
a strongly continuous semigroup lA{t) given by the Trotter product formula 



U{t)x ~ lim 



->C50 



X, xeX, (2.4) 



with uniform convergence for t in compact intervals. 

It follows from the previous result that if T{t) and S{t) are both positive 
then also ZY(t) is positive. It is also clear from (|2.4p that if ||T(t)|| < expjcjii} and 
||5(t)|| < exp{uj2t} then we have ||^^(t)|| < exp{(wi +uj2)t}- 

Our main result is the following. 

Theorem 2.3. The operator A + B generates a positive strongly continuous quasi- 
contractive semigroup {T{t)}t>o of bounded linear operators on X. 

Proof. Our aim is to apply first Theorem l2.1l to show that the operator A generates 
a positive contractive semigroup. To this end, for every u G Dom(^) we define 
0u e X* by 

ui{s) > 0, U2{s) > 0, 

M^) = { ):;:{ :; "^J^Jj^' . (2.5) 

^"v ; (-n iA ,f ui(s) < 0, U2{s) > 0, ^ ' 

ui{s) < 0, U2{s) < 
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Then 

and clearly 



ll0u|U* < 1, 



(u,0u)Ar= f (ui(s)x„+(s) + "2(s)x,4(s)) ds=||u+|U, 



where Xu'^ Xu+ stand for the characteristic function of the support of u'^ and 
U2 , respectively. Making use of (|2.5p we obtain 

{Au,(t>^)x = - I (7i(s)wi(s)),X„+(s)ds - / (72(5)^2(5))^ x„+(s)ds 
Jo Jo 

< - (7i(m)ui(m))x„+(m) - (72(to)u2(to)) X„+ (n^) 

< 0. (2.6) 

Hence the operator A is dispersive. A is also densely defined. It remains to show 
the range condition, i.e. Rg(AX — A) = X. We observe that the equation 

(AJ - ^) u = h (2.7) 

for {hi, /i2) = h E X and A > sufficiently large has a unique solution {ui, U2) = 
u e Dom(^), given by 



ui{s) 

U2{s) 



I - /o ^ M Jo llo ^ d2 
/o ^ dyl /; exp 



72(2) 



^dy' 

72(1/) 



(2.8) 



The fact that u defined by 
hi{s 



is an element of Dom(^) follows from 



< 



7»(s) 
hi{s) 



A 



l^is) Jo 
M' 



exp 



dz 



\htiy)\ 
ii{y) 



dy 



'A: 



1,2, 



for A large enough for some < 00, i — 1,2, that is u e VK^'^(0, to) x M^^'^(0, m), 
and the range condition is satisfied. Theorem 12.11 gives that ^ is a generator of a 
positive contractive semigroup. Since B is bounded and 



B 



|S||X> 



it follows (see e.g. [2 Theorem 1.11, C-II]) that B generates a positive quasicon- 
tractive semigroup and the proof is completed on the grounds of Theorem 12.21 

□ 

Remark 2.4 We rewrite the operator B as 

B = Bi + B2, 

where 

/-(m + ci)^.i 
V -C2W2 



on X. 
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and 

V ci-ui J 

Then we have 

||exp{tS2}|| < cxp{t||S2||} < cxp{t{mB + C)}, 

and 

II cxp{t;Bi}|| < cxp < — t ■ min < inf {/i(s) 4- ci(s)}, inf {c2(s)} 

I I se[0,m] s£[0,m] 



Therefore if 



+ C < min < inf {^{s) + Ci{s)}, inf {02(5)} 

se[0,m] se[0,m] 



then the growth bound loq of the semigroup is negative, hence it is uniformly 
exponentially stable (see e.g. [Sj), i.e. the population dies out. 



3. Asymptotic behavior 

Model (|l.ip - (|1.3[) is linear hence one naturally expects that solutions cither 
grow or decay exponentially in time (unless they are at an equilibrium). In this 
scenario certain stability properties of solutions can be efhciently investigated. For 
example for a simple age-structured population one can often show that solutions 
behave asymptotically as n{a,t) w Ce''*7i*(a), a property called asynchronous ex- 
ponential growth, where r is the Malthusian parameter and n^{a) is the stable 
age-profile. We recall (see e.g. [8]) that in the framework of linear semigroup the- 
ory a strongly continuous semigroup S — {<5(i)}t>o on a Banach space 3^ with 
generator O and spectral bound 

5(0) =sup{Rc(A) : A 6(7(0)} 

is said to exhibit balanced exponential growth if there exists a projection 11 on y 
such that 

lim ||e-^(°)*5(0 -nil =0. (3.1) 

The semigroup S — {S{t)}^yQ is said to exhibit asynchronous exponential growth if 
it exhibits balanced exponential growth with a rank one projection 11. Wc further 
refer the reader to the monographs [21 El Ej for basic definitions and notions used 
throughout this section. Note that balanced exponential growth essentially requires 
that the spectral bound s{0) is a dominant eigenvalue, and that the semigroup S 
is essentially compact, i.e. uiessiS) < s{0), where uiessiS) stands for the essential 
growth bound of the semigroup S (see e.g. [S]). 

Our aim in this section is to carry out a spectral analysis of the governing 
linear semigroup T{t). First we show that solutions of model (|l.ip - (|1.3p exhibit 
balanced exponential growth (with a finite dimensional projector n). In other 
words, model p.ip - p.3p admits a finite dimensional global attractor. Then we show 
that under some further biologically relevant conditions on the model parameters 
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the governing semigroup is also irreducible (see below for the definition) therefore 
establishing that solutions exhibit asynchronous exponential growth, 

Proposition 3.1. The spectrum of A + B can contain only isolated eigenvalues of 
finite algebraic multiplicity. 

Proof. Our aim is to show that the resolvent operator i?(A, A + B) is compact. 
Since B is bounded, it is enough to show that R{X, A) is compact. Since W^'^{0, m) 
is compactly embedded in L^(0,to) by the Rellich-Kondrachov theorem [TJ Theo- 
rem 6.3, Part I], the claim follows on the ground of [HI Proposition II. 4. 25]. □ 
Proposition 13.11 implies that the essential spectrum of A + B is empty. 
This would imply that the governing semigroup T{t) is essentially compact, i.e. 
uJess{T) < s{A + B) holds, if we could estabhsh that the spectral mapping theo- 
rem for the semigroup and its generator holds true and show for example that 
the point spectrum ap {A + B) is not empty. We would like to point out that in 
general the eigenvalue equation 

{A + B- AX)u = 

cannot be solved explicitly due to the transfer terms ci, C2 and due to the very 
general recruitment term. Hence one needs to find an indirect proof to show that 
for the spectral bound s{A -t- S) > — oo holds, indeed. In the remarkable paper [TB] , 
Heijmans used the Krcin-Rutman theorem to give such a proof for a linear size- 
structured model. Here we present a different approach that uses a lower bound 
of the birth process by an operator of rank one. 

Theorem 3.2. The generator A + B has a non-empty point spectrum. 
Proof. We rewrite the Cauchy problem (|2.ip as follows 

^u= + +63 + ^^4) u, u(0) = Uo, (3.2) 
where Bz + Bi = B2 and 

^^^^(^r/?(-,yK(?/)dj/^ on A- 

and 

\ClUlJ 

Fix a nontrivial separable kernel /3* that satisfies 

0</3*(s,y)=/3i(s)/32(y) </3(s,y) 

and denote by B^ the corresponding operator with birth process defined by /?* 
(such a P* exists, since supp/3 is not emapty). The eigenvalue equation 

{A + Bi+B*3- \I)u = 
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admits a non-trivial solution vector u 7^ if and only if A G C satisfies the following 
characteristic equation 

1 = A'(A)^ r^s) r^exp/- r^^fKr)^jir)+y,ir)^ ^^^^^ 



JO Jo 71 (y) I Jy 71 W J 

(3.3) 

It is easily shown that equation (|3.3p admits a unique real solution A, , which is in 
fact a strictly dominant eigenvalue of ^ + Bi + . Bi is linear and bounded and it 
follows from Theorem 1.11 C-II] that it generates a positive semigroup. Then it 
follows from Theorem 12.21 that A + Bi + generates a positive semigroup, since 
B3 is a positive operator. Clearly, B3 — S3 is a positive operator and we obtain by 
Corollary VI. 1.11 in jg] 

-oo<s{A + Bi + B^)<s{A + Bi+B3). (3.4) 

Finally, since B4 is positive and ^ -I- 4- ;B3 is a generator of a positive semigroup, 
we have by the same corollary 

- 00 < siA + Bi + B3) < s{A + Bi+ B3 + Bi) = s{A + B), (3.5) 

and the proof of Theorem 13.21 is complete. □ 



Remark 3.3 We note that K{0) is the net reproduction rate of the population in 
the active phase when the outflow term ci is interpreted as extra mortality (see e.g. 
[2]) and when the birth rate /3 is separable. Hence if K{0) > 1 then this dominant 
eigenvalue A, is positive, while K(0) < 1 implies that A* is negative, finally A* = 
if and only if K{0) = 1, that is, the net reproduction rate of individuals in the 
active phase equals one. 

Remark 3.4 While the spectral bound s{A + B) cannot be obtained explicitly, it 
can be determined approximately by bounding the kernel P (either from above or 
below) by finite sums of separable kernels 

n 
k=l 

Then, instead of (|3.3|) , we obtain the characteristic equation as the determinant of 
an ri by ri matrix and the leading eigenvalue can be determined at least numerically 
in case of concrete model ingredients. 

Lemma 3.5. The semigroup {T{t)}t>o generated by the operator A + Bi +B3 + B4 
is eventually compact. 

Proof. In the view of the Frechet-Kolmogorov compactness criterion in (see 
e.g. [22I Chapter X.l]) we conclude from 

|fi3u(s)-e3u(s')l < / \{f3{s,y)-(3{s',y))\\u,{y)\dy 
Jo 

< ms,-)-(3{s\-)\\^\\u\\i 
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and the continuity of /3 (see assumption (|1.4p ') that the operator S3 is compact. 
Hence it suffices to consider the operator A + Bi + B4. We define the function 7 
as follows 

7(5) = min{7i(.s),72(s)} for s G [0, m], 
then 7 is clearly bounded and strictly positive because of our assumptions on 71 , 72 
in (|1.4p . It follows that i is integrable and 

Tis)=r^dr (3.6) 
Jo lir) 

is greater or equal than the time needed for an individual to grow from size to 
size s irrespective of the transit rates ci and C2 . This implies that the semigroup 
7o(t) generated by the operator A + Bi + Bi is nilpotent, in particular for any 
u e we have 

To(t)u = for t> Tim), 
and the claim follows. □ 

Theorem 3.6. The semigroup T(t) generated by the operator A+B exhibits balanced 
exponential growth. 

Proof. We have already established the existence of a spectral gap for the generator 
A + B, hence it only remains to show that the semigroup generated by ,4 + S is 
eventually norm continuous. Then the Spectral Mapping Theorem holds true and 
it follows (see e.g. Corollary VI. 1.13 in [^) that the boundary spectrum 

(T+{A + B) = g{A + B)r\ {s{A + B) + iM) 

equals s{A-\-B) which is a pole of the resolvent with finite algebraic (hence geomet- 
ric) multiplicity by Proposition 13.11 The eventual compactness of the semigroup 
Tit) established in Lemma 13.51 implies eventual norm continuity, hence the proof 
of Theorem 13. 61 is complete. □ 

Remark 3.7 Here we would like to briefly present another idea to show that the 
semigroup generated hy A + Bi + Bz + B^ is essentially compact. This approach 
works even in the case when individual size may be arbitrary large, i.e. for models 
when s S [0, 00), when one cannot establish eventual compactness of the semigroup 
and the Spectral Mapping Theorem. Since the operator B71, is compact (and positive 
as well) and + Si + S4 is a generator of a positive semigroup we have 

OJess {A + Bi+Bz+Bi)= LJess {A + Bi + B4) < M-^ + ^1 + B4) ^ s{A + Bi + B4) . 

Therefore one only really needs to show that 

s{A + Bi + B4) < siA + B1 + B3 + B4). 

But this follows from Theorem 13.21 and by noting that the eigenvalue equation 

{A + B1+B4- AJ)u = 0, u(0) = 

does not admit a non-trivial solution for any A e C. 

We conclude the section with the following result. 



A size-structured two-phase population model 



11 



Theorem 3.8. Assume that there exists an Eq > such that for all < e < Sq 

I I l3{s,y)dyds>Q (3.7) 

and that the transition rates satisfy 

inf suppci = 0, and sup suppc2 = m. (3-8) 

Then the semigroup T{t) generated by A + B exhibits asynchronous exponential 
growth. 

Proof. It only remains to show that under conditions (|3.7p . p.Sp the semigroup 
T(t) is irreducible, i.e. for every ^ u G and 7^ u* G there exists a to 
such that 

{T{to)u,u*)x>0. (3.9) 
Let u ~ (^1,1*2) and tti and tt2 be the projections onto the first and second 
coordinates, respectively. First assume ui{s) > for s in some interval [sqi 
Since 71 > 0, there exist < t* such that 

supp TTi (T {t)u) n supp /3(s, ■ ) ^ 

for every t^ < t < t* and every s G (0, e]. By assumption (|3.7p . 7ri(T(t)u)(s) > 
for t^, < t < t* and s G (0, e] and eventually in a neighborhood of every s G (0, m]. 
If ui{s) = for all s, then 1*2(5) > for s in some interval [sq , Sq]. Since 72 > 0, 
there exist t* < t* such that 

supp7r2 (T(t)u) n suppc2 7^ 

for every t^ < t < t* and hence by assumption (|3.8p tti (T(t)u) (s) > for s G 
suppc2. By a similar argument as above, tti {T{t)u) (s) > for s G (0,e] and t 
sufficiently large. Finally, assumption p.Sp also guarantees that 7r2 (T(t)u) (s) > 
for s G (0, e] and t sufficiently large. 

□ 

4. Concluding remarks 

Quiescent phases have been introduced in a variety of biological models, see 
e.g. [21 [71 [131 mi and the references therein. They arise in biological populations 
that show distinguished classes of behaviors. For example, it is well known that 
cells in a tissue can either migrate or divide, but not both. To the best of our 
knowledge, a size-structured models with a quiescent phase has not been studied 
yet. 

Asynchronous exponential growth is often observed in linear models of pop- 
ulation dynamics. It represents the fact that although a population is still expo- 
nentially growing (or decaying), it's age- or size distribution is "at equilibrium". 
Several authors have proved the property of asynchronous exponential growth for 
age-structured populations |31 [7], with a possibly delayed birth process [TS]. In 
this paper, using the techniques established for semigroups of positive operators, 
we have proved a corresponding result for a size-structured population with a 
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quiescent conipartment. It should be stated here, that the conditions that distin- 
guish asynchronous exponential growth from merely balanced exponential growth, 
namely (|3.7|) and (|3.8p are natural in the sense that they can be achieved by an 
appropriate choice of the size space. If, for example, the possibility to return to the 
active class would cease at some mi < m, i.e. 02(3) ~ for all s > rrii, it would 
be possible to reduce the size space in the second coordinate to [0, mi], without 
affecting the asymptotic behavior. 

It was shown recently by Hadcler and Thicme |15j for finite-dimensional mod- 
els that coupling to a quiescent phase can shift spectral bound of a matrix (i.e. the 
intrinsic growth rate) in both directions, depending on the transition rates be- 
tween active and quiescent phase. It can be seen easily that the same holds for our 
infinite-dimensional setting. The transition to the quiescent phase can either cir- 
cumvent a region of high mortality or high reproductive activity in the size space 
of the active population. 

We would like to note that the next natural step will be to incorporate 
some interaction variables into our model. Interaction between individuals in the 
population may be induced for example by scramble competition for available 
resources. It also seems reasonable for example to require that the transition rates 
ci and C2 depend on the standing population sizes in the active and resting phases, 
respectively. Then our model becomes a nonlinear one, and the question arises what 
types of nonlinearities will preserve the asymptotic properties of solutions discussed 
in the previous section. We note that results in the literature for nonlinear models 
are rather rare, see for example |12j . 
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